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Abstract 

The Garden of Eden theorems are well known theorems established by Moore and Myhill in the early sixties 
connecting injectivity and surjectivity for the global function of cellular automata in the (Euclidean) plane. 
In this paper, it is shown that the properties established by Moore and Myhill are no more true for cellular 
automata in the hyperbolic plane. 
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1 Introduction 

Cellular automata are a well known model of computation, studied from already fifty 
years. This model is a very powerful tool for simulation of various phenomenas, from 
gas flows up to the circulation of crowds of pedestrians in railway stations and other 
physical, chemical or sociological problems of everyday life. At the same time, they 
are very simple dynamical systems which raise a lot of difficult theoretical questions 
which are far from being solved. Cellular automata are also a model of computations 
to which the theory of computability applies. There are a lot of results about the 
complexity of computations with cellular automata. There are also a lot of results 
on the Turing complete power of computation of this model. 

In the early fifties, Moore and Myhill established an interesting property of 
the global function of cellular automata in the plane, in fact, the Euclidean plane. 
Remember that a cellular automaton A in the plane is defined by the definition of a 
neighbourhood, here the Moore neighbourhood, and a local transition function. For 
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any cell c whose address is given by a couple of integers, (c.x, c.y), c.x : c.y E 'Z, the 
Moore neighbourhood is the set N c of cells d for which | d.x — c.x |, | d.y — c.y \ < 1. 
The local transition function /a is a mapping from Q\ into Qa-, where Qa is the 
finite set of states of A, which rules the change of state in the cell at each tic of the 
clock: 

rf A (c,t+ 1) = f A ({r] A (d,t)\d G N c }), 

where rjA(c,t) is the state of the automaton at the cell c, at time t. 

The global transition function Ga , is defined is a mapping from Q ^ into itself 
defined as follows. If £ E Qa is a configuration, Ga(€) is the configuration at the 
next time and it is defined by: GG(£)( C ) = /a({£(gO} I d E N c }). 

The theorems established by Moore and Myhill, see [20,21] state that for a 
cellular automaton A in the plane with Moore neighbourhood, Ga is surjective if 
and only if it is injective on the set of finite configurations. This property was used 
by Jarkko Kari to establish that the injectivity, the surjectivity and the reversibility 
of cellular automata in the plane are undecidable, see [4]. 

In the nineties and very recently, a few works appeared on the question whether 
the question of the existence of Garden of Eden is connected with the growth func¬ 
tion of the group associated to a Cayley graphs. This function of n gives the number 
of elements of the group which can be represented by a product of at most n gen¬ 
erators of the group and their inverses. In [6,2,1], the question is considered in 
different settings of Cayley groups. In [6], a counter-example is given when the 
group has an exponential growth: then, there is an injective global function which 
is not surjective and there is a surjective global function which is not injective. 

However, the techniques of Cayley graphs do not apply to infinitely many tes¬ 
sellations of the hyperbolic plane. They apply only when the generating polygon 
is of the form {2p, 2g}, which means that its number of sides is even and that the 
number of polygons around a vertex in the tiling is also even, see [3]. In the other 
cases, especially for polygons with an odd number of sides and with an odd num¬ 
ber of copies around a vertex, these techniques do not apply. However, for all the 
mentioned tessellations, whatever the parities, the technique indicated in [7] always 
apply: there, tools are developed which allow to locate the cells in the hyperbolic 
plane, facilitating the implementation of cellular automata in these contexts. 

Cellular automata in the hyperbolic plane appeared for the first time in [16,17]. 
After the publication of [7], this study received a new impulse, see for instance [11] 
for references. In [10], I have shown that Hedlund’s characterization of cellular 
automata in the plane can be transported to cellular automata in the hyperbolic 
plane at the price of an additional property on the cellular automaton. I also have 
noticed that the proof of Moore and Myhill in the Euclidean case does not pass 
to the hyperbolic plane. In the light of the mentioned paper in the case of Cayley 
graphs, it seems plausible that the answer is that there are also counter-examples. 
In fact, as shown in Section 3, there are examples of an injective cellular automaton 
which is not surjective and of a surjective cellular automaton which is not injective, 
and the local transition function which we give in these examples are much more 
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simple than in [6]. 

In Section 2, we remind a few notions about cellular automata in the hyperbolic 
plane, also see [11]. 

2 Cellular automata in the hyperbolic plane 

It is not needed to be very familiar with hyperbolic geometry to have a good rep¬ 
resentation of what happens in the hyperbolic plane. Fortunately, there are good 
models for that. In this paper, we make use of Poincare’s disc model. The hyper¬ 
bolic plane is the set of points inside a fixed open disc U of the plane and the points 
of dU, the border of £/, are called the points at infinity. In this model, lines are the 
trace in U of diametral lines or of circles which are orthogonal to dU : see Fig. 1, 
below. The interest of this model is that the angles between lines in the model are 
the true angles in the hyperbolic plane. 



Fig. 1. The Poincare’s disc model. Remark the two lines p and q. They both pass through the point A 
and they are both parallel to the line £. The line m passes through A and does not cut the line £, even at 
infinity. 


Cellular automata live on regular grids. There are infinitely many ones in the 
hyperbolic plane. Here, we choose the simplest ones in some sense: the pentagrid 
and the ternary heptagrid, see [11] for a detailed study of these grids. The pentagrid 
is constructed on the replication of a copy of the regular pentagon with right angles 

in its sides and, recursively, of the images in their sides. The ternary heptagrid 

2 7T 

is constructed in a similar way starting from the regular heptagon with — as its 
angles. 

A view of these grids are given in Fig. 2 in a way which indicates how it is 
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possible to locate cells and which plays a key role in the proofs of the theorem of 
Section 3. 

The important property indicated by Fig. 2 is that in the case of the hyperbolic 
plane, the penta- and the heptagrid are generated by a tree. The tree structure 
is underlined by the parts of the tiling which are detached and placed around the 
central tile. In the case of the pentagrid, five such regions, each one spanned by 
the same tree, are placed around the central tile, in a rotation symmetric way. In 
the case of the ternary heptagrid, we have seven regions. A remarkable property is 
that the generating tree is the same for the pentagrid and for the ternary heptagrid, 
see [11]. This tree is called Fibonacci tree as the number of its cells on a level k 
is / 2 /c+i, where {fk}keiN is the Fibonacci sequence where /o = fi = 1, see [16]. 
From this, coordinates can be computed to locate the cells of a cellular automaton, 
see [7,11]. 



Fig. 2. On the left: the pentagrid, on the right: the underlying tree which spans the tiling. 


As a consequence, except the cell which is placed at the central tile, we call it 
the central cell, each cell of the cellular automaton, has a father: its father as a 
node of the tree in the region it falls in. 

This is an important point which will be used in Section 3. 

The existence of a father for all cells, except the central one, plays the role of a 
direction in the hyperbolic plane, in the same way as the four traditional directions 
play a key role in the Euclidean plane. 

In most formal presentations of cellular automata in the Euclidean plane — 
people usually say CA in the plane — the set of cells is identified with 'Z 2. This 
identification is so evident that it requires some effort to realize that it connects two 
different things and how it performs the connection. 

In the Euclidean case, the above identification consists in three steps. First, we 
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fix a cell whose coordinates will be (0,0), by definition. Then, we fix the directions, 
North , South , East and West and this requires to fix two cells: the one with 
coordinates (1,0) and the one with coordinates (0,1): the second choice amounts 
to define what means clockwise. These three choices are arbitrary as the Euclidean 
plane has no privileged point and as it has no intrinsic orientation. 

Now, the hyperbolic plane has no privileged point and it is intrinsically non- 
oriented too. Basically, we do things in a very similar way with what is performed 
in the Euclidean case. First we fix a cell which will be the central one. This fixes 
the sectors around the central cell: each neighbour of the central cell is identified 
with a root of the Fibonacci tree which spans the considered sector. Next, fixing 
in a sector which cell will be the leftmost son of the root allows us to fix what 
means clockwise in this context, as the leftmost son of another sector is defined by 
the rotation which maps this sector to the first chosen one. This correspondence 
between the two processes allows us to say that defining the father of all the cells, 
except the central one, amounts to define a direction in the hyperbolic plane. 

3 Gardens of Eden in the hyperbolic plane 

Now, we turn to the examples which we announced in the introduction. 

Theorem 3.1 (Kari-Margenstern) — There is a cellular automaton A on the 
ternary heptagrid, or on the pentagrid, such that Ga is injective but Ga is not 
surjective. There is also a cellular automaton B on the ternary heptagrid, or on the 
pentagrid, such that Gb is surjective but Gb is not injective. 

Proof. First, consider the case of A. 

We assume that the ternary heptagrid, or the pentagrid, has coordinates based 
on a central cell xq and the required number of sectors around it, each sector being 
spanned by a Fibonacci tree, as mentioned in Section 2. 

Each cell x, with x ^ xo, has a father which we denote by f(x). For A , we 
consider that there are two states, 0 and 1 and that the transition function is 
defined by the following relations: 

r] A (x 0 ,t+l) = rj A (x 0 ,t), 

rj A (x,t+l) = rj A (f(x),t), when x + x 0 , 

where rjA(y,t) is the state of the cell y at time t under A. 

It is not difficult to see that Ga is injective. Indeed, if c\ and C 2 are two 
configurations with c\ ^ C 2 , there is a cell x such that c\(x) ^ C 2 (x). If we consider 
a son y of x, we have that Ga(ci)(v) = c\{x) ^ C 2 (x) = Ga(c 2 )(y)- Now, Ga cannot 
be surjective as it cannot reach configurations in which two sons of the same node 
have different states. 

Let us turn to the construction of B. We can give two variants of this construc¬ 
tion. 

In the first variant, we assume that we have a function ac which, to each node, 
associates its leftmost son. If x ^ xq, where xq again denotes the central cell, cr^(x) 
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is known by the cell from f(x). For xo, as the sectors spanned by a Fibonacci tree 
are numbered, we define cr^(xo) as the root of the tree which received the smallest 
number. 

Now, we define B as follows: 

VB(x,t+ 1 ) = xoT(rj B (x,t),riB(<Te(x),t)), 

where rjB(x,t) is the state of the cell x at time t under B. 

It is not very difficult to see that Gb is not injective. If we define Co by assigning 
the state 0 to all cells and c\ by assigning the state 1 to each cell, it is not difficult 
to see that G b (cq) = G B (c i) = cq. 

Now, let us check that Gb is surjective. Indeed, fix a configuration c\ and we 
have to define a configuration Co such that G b (cq) — c\. 

Consider xq. Define co(xo) = 0. Then, applying the definition, we have that 
Ci(x 0 ) = xor(0, co(cr^(xo))) = co(cr^(xo)). And so, this defines Co at &i(x o). Define 
c 0 (x) = 0 for all the other sons of xq than a^(x o). 

By induction, assume that we have defined the level n+1 and that the surjectivity 
holds for all cells up to the level n, this level being included. From what we have 
seen, this is the case for n = 0. 

On the level n+2, define all white nodes y by co(y) = 0. Now, consider a node x 
of the level n+1. As co(x) = a is defined, we have, by definition, Gb(cq)(x) = 
xor(a, co(cr^(x))). This always defines Co at cr^(x). Indeed, co(cr^(x)) = ci(x) if 
a = 0 and co(cr^(x)) = 1—ci(x) if a = 1. And so, considering all nodes x of the 
level n+1, this defines Co for all black nodes of the level n+2. Moreover, now 
Gb(co)(x) = c\{x) for all cells x of the level n+1 too And the definition of Co at the 
level n+2 is complete. 

And so, by induction, we proved that G B is surjective. 

The second variant requires to know the sons of a cell. This is easy to define 
from /(x) for any cell x with x ^ xq. For xo, we consider that all its neighbours 
are its sons. Now, we define as follows: 

r] B (x,t+ 1) = mm{x.or(r]B(x,t),r]B(y,t)) \y £ S x }, 

where S x is the set of the sons of x. 

The argument is the same as in the first variant. 

Note that, in the proof of the surjectivity, we can easily see that B cannot be 
injective, as long as the state of many cells can be fixed arbitrarily. □ 


4 The case of rotation invariant cellular automata 

It was proved in [15] that an analog of Hedlund’s for cellular automata hold for the 
hyperbolic plane provided that an additional property is satisfied by the automaton, 
namely that the set of its rules is rotation invariant. 

Intuitively, this means that if the neighbourhood of a cell is changed by a rotation 
of the neighbourhood around the cell, then the new state of the cell is the same as 
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what it was defined before the change. This informal definition can be made more 
precise, by fixing a way to number the neighbours of a cell or, which is equivalent, 
to number the sides of the polygon which supports the automaton. Let us say that 
for the central cell, side 1 is fixed once and for all. The other sides are numbered 
increasingly while counter-clockwise turning around the cell starting from side 1 . 
For the other cells, fix side 1 to be the side shared by the cell with its father and, 
similarly number the other sides by counter-clockwise turning around the cell. Then, 
a rule of the automaton can be displayed in the following format, see [12,19,18] : 

m,vi,-,Va -> vl, 

where 770 is the current state of the cell, 77 ^, with i E {l..<a}, where a = 5 or a = 7, 
is the current state of neighbour i and t?q is the new state of the cell. 

We say that a cellular automaton on the pentagrid or on the heptagrid is rotation 
invariant if and only if for each rule as above, the rule 

VOi ^tt(I) j •••? r hr(a) > 

also belongs to the set of rules for any circular permutation n on {1 ..a}. 

In what follows, s(x,t) denotes the state at x and at time t. We have the 
following two properties : 

Theorem 4.1 (Margenstern-Kari) There is a rotation invariant cellular automa¬ 
ton on the pentagrid or on the heptagrid which is surjective but not injective. 

Theorem 4.2 (Margenstern-Kari) There is a rotation invariant cellular automa¬ 
ton on the pentagrid or on the heptagrid which is surjective but not injective even 
on finite configurations. 

Proof of Theorem 4.1. It will be enough that s(x,t) takes its values in {0,1}. 
We define the transition rules in such a way that: 

, (O if #{s(y,t) = 0 } is even 

s(x, t + 1 ) = s , 

[l if #{s( 77 , t) = 1 } is even 

where, in this formula, y denotes a neighbour of x which is not x and t))} 

is the number of t/’s around x such that s(y,t) satisfy P. This is clearly a rotation 
invariant function. 

Consider a configuration c and again, denote the central cell by xq. It is not 
difficult to see that co(xo) can be 0 or 1 indifferently. Once it is fixed, we can fix 
the level 1 around xo, setting two cells at the value c(x 0 ) and the others at 1 —c(xo). 
Then, by induction, on each level, we have that the number of neighbours to be 
fixed is 3 or 2. In fact, the number is 4 for the first tile which we consider at level 1. 
For the others it is 3, the last one being excepted, for which it is 2. Now, having 2 
or 3 cells at our disposal is enough to fix the parity of the required value to be odd. 
If the number of the required value is odd it is enough to give this value to one new 
cell and the other value to the others. If it is even, we have the choice: even to give 
the required value to an even number of cells, 0 or 2 , we have the choice, and the 
other value to the other cells. 

This proves that the global function is surjective. Moreover, the just above 
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argument shows that it is not injective: whatever the number of cells, we have a 
choice. Also, we already have a choice for co(xq). □ 

Proof of Theorem 4.2. The states are again 0 or 1 and we denote by N(x,t, 1) 
the number of neighbours y of x : with y ^ x for which s(y,t) = 1 , where, again, 
s(y 1 t) is the state at y and at time t. We define the transition function in such a 
way that: 

(0 if N(x,t,l) e {0,1,4,7} 
s(x, t + 1) — \ , 

[1 if N(x,t, 1) e { 2 , 3, 5, 6 } 

It is again clear that this transition function can be defined by rotation invariant 
rules. 

The surjectivity comes from the fact that if we have two cells at our disposal 
among the neighbours of a cell x, this is enough to fix the value at x according to 
this rule. This can be checked by Table 1, below. 

In the table, N*(x, £, 1) is the number of l’s on the already fixed neighbours of x. 
Each entry tells us how many cells to put to 1. If the entry says 0, this means that 
both free cells are put to 0 . 

Note that if we have more free cells at our disposal, we can use the same table 
to fix the values: two cells a and b among the free ones are fixed according to the 
table and the others are fixed to 0 or to 1 in order to obtain the value 7V*(x, £, 1 ) of 
the table for the cells which are distinct from a and b. 

Table 1 The values of N* (x,t, 1) for the example of Theorem 4.2. 


1 ) 

0 

1 

2 

3 

4 

5 

if 0 needed 

0 

0 

2 

2 

0 

2 

if 1 needed 

2 

2 

0 

0 

2 

0 


Now, the configuration co where every cell has the state 0 is transformed into Co 
and the configuration c\ where every cell has the state 1 is also transformed into co- 
Accordingly, the global function is not injective. Now, it is also not injective on 
finite configurations. 

Indeed, let us fix a cell xq. The configuration c(x o) defined by 1 at xo and 0 
everywhere else is also transformed into cq. Similarly, if we have scattered l’s at a 
distance at least 4 from each other, then such a configuration is also transformed 
into co- Indeed, the requirement on the distance entails that the neighbourhoods of 
the l’s are disjoint. And so, we can take only finitely many of them if required. □ 

5 Conclusion 

We have proved that injectivity and surjectivity are independent in the case of 
cellular automata in the hyperbolic plane. In the case of rotation invariant cellular 
automata, the problem of finding an example of an injective cellular automaton 
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which would not be surjective is still open. 

From these results, it appears that a new argument has to be found to prove 
the undecidability of the surjectivity or the reversibility of cellular automata in the 
hyperbolic plane. For the injectivity, the problem can be proved undecidable by 
transferring to the hyperbolic plane the way opened by Jarkko Kari, although it is 
not all that easy, see [13]. It seems reasonable to conjecture that the surjectivity and 
the reversibility of the global function of a cellular automaton are undecidable. An 
argument, in favour of this conjecture is that the tiling problem for the hyperbolic 
plane is also undecidable, as it was proved by Maurice Margenstern and by Jarkko 
Kari, independently and by very different methods, see [5,8,14]. 
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